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ABSTRACT 

In this paper, we describe a novel experiment for the accurate estimation of pulsar 
dispersion measures using the Giant Metre-wave Radio Telescope. This experiment 
was carried out for a sample of twelve pulsars, over a period of more than one year 
(January 2001 to May 2002) with observations about once every fortnight. At each 
epoch, the pulsar DMs were obtained from simultaneous dual frequency observations, 
without requiring any absolute timing information. The DM estimates were obtained 
from both the single pulse data streams and from the average profiles. The accuracy 
of the DM estimates at each epoch is ~ 1 part in 10 4 or better, making the data set 
useful for many different kinds of studies. 

The time series of DM shows significant variations on time scales of weeks to 
months for most of the pulsars. A comparison of the mean DM values from these 
data show significant deviations from catalog values (as well as from other estimates 
in literature) for some of the pulsars, with PSR B1642— 03 showing the most notable 
changes. From our analysis results it appears that constancy of pulsar DMs (at the 
level of 1 in 10 3 or better) can not be taken for granted. For PSR B2217+47, we see 
evidence for a large-scale DM gradient over a one year period, which is modeled as 
being due to a blob of enhanced electron density sampled by the line of sight. For 
some pulsars, including pulsars with fairly simple profiles like PSR B 1642— 03, we find 
evidence for small changes in DM values for different frequency pairs of measurement, 
a result that needs to be investigated in detail. Another interesting result is that we 
find significant differences in DM values obtained from average profiles and single pulse 
data. 

Key words: miscellaneous - methods:data analysis - pulsars: general - Hn regions. 



1 INTRODUCTION 

The radio signals from a pulsar suffer dispersion as they 
travel through the ionized component of the inter-stellar 
medium (ISM), resulting in a frequency dependent arrival 
time of the pulses. The effect is quantified by the pulsar's 
dispersion measure (DM), defined as the integral of the elec- 
tron column density along the line of sight, 



DM = 



riedl pc/cm 



(1) 



The delay between the pulse arrival time at two frequencies, 
At, can then be expressed as 



At = K 



where 



1 

7F 



/I 



DM 



K = 



2mnc 2.410331 x 10" 4 



MHz 2 cm 3 s/pc 



(2) 



(3) 



Here At is in units of second for fi and ji in MHz and DM 
in the traditional units o f pc/cm 3 . The precis e value of the 
constant K is as given in lBacker et al. I lll993F ). 

The DM of a pulsar is a basic parameter, and its value 
needs to be known with sufficient accuracy for proper dis- 
persion correction to be carried out on the received sig- 
nal. Further, accurate estimates of DM can be used to 
prob e the pulsar emission geometry (e.g. iKardashev et al~l 
1982). Estimates of DM obtained from different values of 
/i and fi in Equation [3 have been used to check the valid- 
ity of the cold plasma dispe rsion relation for the ISM (e.g. 
IPhillios fc Wolszczan lll992l and references therein). In ad- 
dition, small variations in a pulsar's DM are expected due 
to random electron density fluctuations in the ISM, thought 
to be associated with turbulence in the medium. Such vari- 
ations, expected on relatively large tim e-scales of weeks to 
months, have indeed been ob served fe.g. lBacker et al. Il993l : 
IPhillins fc Wolszczan I Il99lh . Pulsar dispersion monitoring 
thus provides a direct method for probing the structure of 
the spectrum of electron density fluctuations. 
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Though first order estimates of the DM can be ob- 
tained by careful measurements of the arrival time de- 
lays in a multi-channel receiver operating at a single wave- 
band (e.g. during the pulsar search and discovery process 
itself) the more accurate estimates needed for the appli- 
cations discussed above require more sophisticated exper- 
iments. Typically, refined pulsar DMs (and their variations 
with epoch) are estimated as part of the analysis of multi- 
epoch multi-frequency timing data from an obser vatory (e.g. 
iBacker et aT"lll99.i : Iphillips fc Wolszczan Nl992T l. An alter- 
nate method is to conduct simultaneous dual frequency ob- 
servations at /i and /2 and estimate the DM fr om a measure 
of the arrival time delay, using EquationCO(e.g . lBartel et al. I 
Il98lt iKardashev et al. lllQsi lHankins fTlflSTT I. The advan- 
tage of this method is that observations at a single epoch 
are self-sufficient for obtaining the DM at that epoch and 
the DM is obtained more directly, rather than as one of the 
parameters in a multi-parameter timing solution. For single 
dish telescopes, this method requires simultaneous opera- 
tion of receivers at more than one wave-band; alternatively, 
different single dish telescopes can be configured at each 
wave-band while simultaneously observing the same pulsar. 

In this paper, we describe a new experiment for accu- 
rate estimation of pulsar DMs, using the Giant Metre-wave 
Radio Telescope (GMRT) in a simultaneous multi-frequency 
pulsar observation mode. Section [2] describes the details of 
the experiment and the observation strategy. SectionElgives 
the details of the data reduction, and describes the technique 
used for estimating DMs from the reduced data. The main 
results and the possibilities for follow-up work are described 
in Section [1] 



2 A NEW EXPERIMENT FOR MEASURING 
DM 

The accuracy of the DM estimate depends on the precision 
to which the the time delay between the pulse profiles at 
two frequencies can be measured. If At rms is the error on 
the measurement of the time delay, then the fractional DM 
error is 

DA'l rr ns Atrms / . \ 

DM = At ■ [ ' 

For a given value of At rms (which is usually limited by the 
S/N of the data at the two frequencies, or sometimes by 
the coarseness of the sampling interval), it is clear that the 
greater the relative time delay between the arrival of signals 
at the two frequencies, At, the more accurate is the DM es- 
timate. This would favour large separations between the two 
observing radio bands. However, if the pulsar profile evolves 
significantly over this range of frequencies, then it can bias 
the measured At, leading to an error in the estimate of the 
DM. This effect favours a smaller separation between the 
two radio wave-bands. Also, according to Equation for a 
given separation between a pair of radio bands fi and /2, 
smaller values of frequencies give a larger value of estimated 
At, and in turn, a better accuracy for final DM estimation. 
The final choice of the two frequency bands of operations 
is then decided by these considerations. Other requirements 
for obtaining accurate DM estimates are (a) high signal to 
noise ratio stable pulse profiles, which are more readily ob- 
served at low radio frequencies (typically in the range 100 



to 1000 MHz) where the pulsar is known to be bright and 
(b) accurate time alignment of the multi-frequency pulse 
profiles. As we now describe, the GMRT, because of some 
unique features, offers a novel way for obtaining accurate 
DM estimates. 

T he GMRT is a multi- element aperture synthesis tele- 
scope l|Swarup et al. lflQ97fl consisting of 30 antennas, dis- 
tributed over a region of 25 km diameter, which can also be 
configured as a "single dish" in the incoherent or coherent 
array mode flOupta et al. l2o'()fT l. Furthermore, it supports a 
"sub-array" mode of operation where different sets of anten- 
nas can be configured completely independently to produce 
more than one single dish. Thus, the same pulsar can be 
observed simultaneously at more than one radio band. 

The GMRT operates at radio frequencies in the range 
150 MHz to 1400 MHz with observing bands available at 
150, 235, 325, 610 and 1400 MHz. The antennas can be 
grouped into several sub-arrays and each sub-array can in- 
dependently be operated at a radio band of interest, thus 
enabling simultaneous multi-frequency observations. Signals 
from different observing frequency bands and antennas are 
eventually down-converted to baseband signals of 16 MHz 
band-width. The signals are subsequently sampled at the 
Nyquist rate and processed through a digital receiver sys- 
tem consisting of a correlator and a pulsar back-end. 

For each antenna, operating at a given frequency band, 
the pulsar back-end receives signals in 256 channels span- 
ning the band-width of 16 MHz, for each of two orthogo- 
nal polarizations. The relative delay - geometrical as well 
as instrumental - between different antenna signals is com- 
pensated to an accuracy of 32 nanosec before they reach 
the pulsar receiver. The corresponding signals from selected 
antennas (say from one sub-array) can be added together 
incoherently by the pulsar receiver. 

For this experiment, the signals from antennas in all 
sub-arrays were added incoherently in the same pulsar re- 
ceiver, to produce a single stream of output data, which was 
recorded at a sampling rate of 0.516 millisecond. Because 
of the dispersive delay between the different radio bands of 
observation, the pulse arrives at different times (and hence, 
at different pulse phases) at each frequency band. This fact 
is utilised to separately extract the streams of single pulses 
at each frequency band, from the single stream of recorded 
data, during the offline analysis. This scheme eliminates 
the need for having separate, but synchronised, pulsar re- 
ceiver chains for each sub-array and also does away with 
any requirement of accuracy of absolute time stamping of 
the recorded data - the data from the different sub-arrays 
is naturally synchronised. Since all known instrumental and 
geometric delays have been corrected for all the sub-arrays, 
the residual arrival time delay between pulses from different 
radio bands of observation is only and entirely due to the 
dispersion delay. This allows the DM to be measured to a 
very high degree of accuracy. 

There is, however, one drawback of the above scheme. In 
order to recover the pulsed signal for the different frequency 
bands during off-line analysis, dispersion delays across the 
256 channels (16 MHz baseband band-width) for each fre- 
quency band are computed and the data are collapsed to 
obtain a time series for each band. In this process, how- 
ever, the data from the other frequency band are wrongly 
de-dispersed and appear as a smeared out signal producing 
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Table 1. Relevant parameters of our selected sample of pulsars. 



Pulsar Catalog DM Period S400 Distance V pm Duration of Frequency combination 

(pc/cm 3 ) (sec) (mjy) (kpc) (km s -1 ) scan (min) of observation (MHz) 



B0329+54 


26.776 





7145 


1650 


1 


43 


145 


33 


227 


243 - 


b 610 


626 


B0818-13 


40.99 


1 


2381 


100 


2 


46 


376 


22 


227 


243 - 


- 325 


341 


B0823+26 


19.4751 





5307 


65 





38 


196 


22 


227 


243 - 


r 325 


341 


B0834+06 


12.8579 


I 


2738 


85 





72 


174 


22 


227 


243 - 


b 325 


341 


B0950+08 


2.9702 





2531 


400 





12 


21 


22 


325 


341 - 


b 610 


626 


B1133+16 


4.8471 


I 


1877 


300 





27 


475 


33 


325 


341 - 


b 610 


626 


B 1642 -03 


35.665 





3877 


300 


2 


90 


660 


11 


325 


341 - 


b 610 


626 


B1642-03 


35.665 





3877 


300 


2 


90 


660 


11 


227 


243 - 


b 325 


341 


B1919+21 


12.4309 


1 


3373 


200 





66 


122 


11 


227 


243 - 


- ,325 


341 


B1929+10 


3.176 





2265 


250 





17 


86 


11 


227 


243 - 


b 325 


341 


B1929+10 


3.176 





2265 


250 





17 


86 


22 


325 


341 - 


b 610 


626 


B2016+28 


14.176 





5579 


320 


1 


10 


12 


11 


227 


243 - 


b 314 


320 


B2016+28 


14.176 





5579 


320 


1 


10 


12 


22 


325 


341 - 


b 610 


626 


B2045-16 


11.51 


I 


9616 


125 





64 


289 


11 


227 


243 - 


b 314 


320 


B2217+47 


43.54 





5385 


135 


2 


45 


375 


22 


325 


341 - 


b 610 


626 
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50 100 150 200 250 300 
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Figure 1. Dispersion curves across the 16 MHz of base-band 
signal for pulsar B2016+28. The upper panel shows (left to right) 
the dispersion curves for the 243 to 227, 325 to 341 and 610 to 626 
MHz bands of observation & the lower panel shows the dispersion 
curves for 243 to 227, 320 to .304 and 610 to 626 MHz bands, 
respectively. The dotted curves on both sides of the continuous 
curves delineate the extent of the 50% width of the average profile. 



excess undesired power in the off-pulse region. In some cases 
this may overlap with the on-pulse signal from the desired 
frequency band, resulting in corruption of the data. Thus, 
in order to obtain undistorted signals, it is essential that 
we choose an observing strategy that avoids such overlaps. 
This requires us to examine the detailed nature of the DM 
delay curve at each frequency band of interest, and to en- 
sure that the curves do not intersect each other within the 
16 MHz of baseband band-width. In Figure U we show an 
example of this. Here, the upper panel (from left to right) 
shows the dispersion curves for pulsar B2016+28 in the fre- 
quency bands 243 to 227, 325 to 341 and 610 to 626 MHz, as 



seen in the base-band signal, after removal of all delays that 
are integer multiple of the pulsar period. It shows that the 
two dispersion curves at frequency bands 243-227 MHz and 
325-341 MHz intersect with one another over certain range 
of channels. Hence, this combination of frequency bands can 
not be used for such observations of this pulsar. By suitably 
changing the value of the local oscillator signals used for the 
down conversion of the radio frequency bands to base-band 
signals, the range as well as the direction of the radio fre- 
quency signals that span the 16 MHz band-width can be 
changed, thus ensuring proper separation of the dispersion 
curves. In this particular case, it has been achieved by mov- 
ing the local oscillator such that the 325 MHz band covers 
320 to 304 MHz (see lower panel of Figure [y. Appropriate 
frequency combinations were found for each pulsar in our 
sample. 

For this experiment, we selected a sample of 12 pulsars 
having sufficiently large fluxes (5400 > 100 mjy), a range 
of DM values (~ 10 — 40 pc/cm 3 ), and sampling different 
directions in the Galaxy. The relevant parameters are sum- 
marised in Tabled where columns 2,3,4,5 and 6 give the 
values of the DM, period, flux at 400 MHz, distance and 
proper motion respectively, as obtained fro m the pulsar cat- 
alog of iTavlor. Manchester fc Lvne I il993h . At every epoch 
of observation, each pulsar from our sample was observed for 
a few thousand pulses (column 7 gives the duration of the 
observing scan) at a pair of frequency bands (given in col- 
umn 8 of Table selected from the available bands of the 
GMRT. The epochs were separated by intervals of about 
two weeks, and the whole experiment was carried out over 
a duration of about one and half years. 



3 DATA REDUCTION AND ESTIMATION OF 
DM 

The recorded data were pre-processed off-line to convert 
from raw time-frequency format to a single pulse time se- 
ries and folded profiles. The pre-processing involved de- 
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Time sample 

Figure 2. Average pulse profiles of the pulsar B1642— 03 observed 
at 610 (solid curve) + 325 (dotted curve) MHz bands combina- 
tion. The upper and lower panels show the pulse profiles before 
and after the alignment respectively. The excess power regions 
near both edges of the profile at 610 MHz are examples of wrongly 
de-dispersed data from the other band. 



dispersion of the data in two frequency bands, folding and 
interference rejection. 

For each pulsar, to recover the pulse trains at the two 
radio bands, the acquired data were de-dispersed within the 
16 MHz band-width of each band by using the catalog DM 
values given in Tabled Where needed, bad data points were 
rejected from the de-dispersed data. For this, after mask- 
ing the data from the on-pulse regions, the running mean 
data from the off-pulse regions was computed and subtracted 
from the original data. Next, off-pulse data points with am- 
plitude greater than the threshold value (typically chosen as 
3 times the off-pulse RMS) were flagged. In addition, data 
were scanned visually, and manual editing of bad data due to 
radio frequency interference was carried out, where needed. 
At the end of the data rejection step, if a large fraction of 
the data around any on-pulse window was found to be bad, 
the entire pulse was flagged. 

The de-dispersed, interference free data trains were 
folded at the Doppler-corrected pulsar periods to obtain the 
average pulse profiles at the two radio frequency bands (see 
Figure |21 for an example). The pulse profile data at each 
observation band were demarcated with three windows — 
two off-pulse and one on-pulse window. The on-pulse win- 
dow contained the properly de-dispersed average pulse pro- 
file, while the off-pulse windows (one on each side of the 
on-pulse) were off-pulse regions which were free of contam- 
ination from the wrongly de-dispersed pulse profile of the 
other frequency band. Data only from these window regions 
were used in the subsequent analysis described below. 

From the reduced data, the dispersion delay between 
the two frequency bands was estimated and, using Equa- 
tion|2| the corresponding DM value was obtained. For these 
calculations, Doppler corrected frequencies /i and /2 (with 



fi > fz) were used, with these frequencies being related to 
the frequencies of observations, /i m and /2m , through 

h = f im \Jjzrji and = f 2m J r^rf ' 13 = > ( 5 ) 

where v ne t is the radial velocity of the observer with respect 
to the pulsar, which is predominantly due to the orbital 
motion of the earth around the Sun. Similarly, the value of 
At in Equation[2]needs to be the measured topocentric delay, 
Atm, corrected to the solar system barycenter, as follows: 

At = Atm X (1 - 0) . (6) 

The total measured time delay, At m , can be expressed as a 
sum of three terms: 

At m = At p + AU + At f , (7) 

where At v is the integral number of pulsar periods delay, 
Ati is the number of time sample bins delay within a pulsar 
period and Atf is the fraction of a time sample bin delay. 
The value of At m can be estimated by two different tech- 
niques: (i) by estimating the delay between the average pulse 
profiles, and (ii) by measuring the mean delay between the 
single pulse data trains. We have carried out the analysis us- 
ing both these methods, and the steps for each are described 
below. 

As the first step, the data were reduced to zero mean 
off-pulse sequences. In the average profile (hereafter AP) 
method, the mean from the off-pulse data windows was es- 
timated and subtracted from the whole pulse profile data. 
In the single pulse (hereafter SP) method, the mean compu- 
tation and baseline subtraction was carried out individually 
for each pulse, while using the same off-pulse windows. 

In the AP method, because of the folding process, the 
value of Atp can not be directly estimated from the folded 
profiles; instead, it was estimated from the knowledge of 
the frequencies for the two bands, the catalog DM value 
and the pulsar period. To estimate Ati, pulse profiles at the 
two frequency bands were cross-correlated, and the integer 
time sample lag at which the cross-correlation peaked was 
taken as value of Ati. The lower frequency pulse profile was 
rotated left circularly by this amount to align it with the 
higher frequency pulse profile (see Figure|3for an example). 

The cross-correlation (hereafter CC) of the pulse pro- 
files at two given frequencies (see Figure |^1 for an example) 
can be given as, 

JV 

CC(kT) = ^f{nT)g(nT-kT). (8) 

71=1 

Here, CC(kT) is the CC for k th bin shift of the pulse profile 
at the lower frequency, N is the number of time sample bins 
within an on-pulse window, and / and g are the pulse pro- 
files at two observation frequencies. In the SP method, the 
two time series were cross-correlated, and the peak of cross- 
correlation function gave the time delay with an accuracy of 
a time sample bin. In this method, the CC could be started 
from zero shift of the lower frequency pulse profile, but to 
reduce unnecessary computations, we started CC computa- 
tions from a shift equivalent to the number of time sam- 
ple bins corresponding to At p . During the cross-correlation 
computations in both the methods, care was taken to en- 
sure that data points from the wrongly de-dispersed signals 
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890 900 910 

Shift (in ms) 

Figure 3. The normalized cross-correlation function (CCF) for 
pulsar B1642-03 observed at 610+325 MHz bands. The contin- 
uous curve shows the CCF for average pulse profiles and dashed 
one correspond to single pulse analysis. 




-20 -10 10 20 -20 -10 10 



Frequency bin Frequency bin 

Figure 4. Normalized CS amplitude (upper panels), and CS 
phase with error bars (lower panels) of average profiles (left side 
panels) and single pulses (right side panels) for pulsar B1642— 03, 
at one epoch observed at 610+325 MHz bands. The straight line 
in the phase plot is the best fit linear gradient. 



were not included in the computations. This was done by 
using data points from the above defined on-pulse and off- 
pulse windows only, and restricting the lag range to values 
which ensured no overlap of these windows with wrongly 
de-dispersed data points. 

The average profile is obtained by folding the time series 
data at the pulsar period. Since individual pulses show sig- 
nificant pulse to pulse jitter in the longitude of occurrence, 
the average profile is usually significantly broader than the 
individual pulses. As a result the CCF obtained in the AP 
analysis is broader in comparison to that from the SP anal- 
ysis (e.g. Figured!. In the AP method, the CCF reflects the 
sum of cross-correlation of all pulses at one radio band with 
all pulses from the other band, while in the SP analysis, the 
CCF is the sum of the CC between corresponding pulses at 
the two radio bands. Therefore, one can expect the DM de- 
lay estimated by the two methods to be different, as we find 
in our results. 

The precision of DM measurement mainly depends on 
the accuracy in estimating the time delay between two pulse 
profiles. The CC as described above gives an accuracy of the 
order of an integral time sample bin. To estimate the delay 
with an accuracy of a fraction of a time sample bin, the 
cross-spectrum (CS) was computed and a linear gradient 
was fitted to the phase of the CS. Let us first consider the 
AP method. If the two pulse profiles are f(t) and g(t), then 
their Fourier transforms (FT) can be written as, 



/(*) 
and 

g(t) 



F{v) 



\F{v)\e 



G(u) = |G(^)|e l (^>+2^2/) 



(9) 



(10) 



where and |G(i/)| are the amplitudes of Fourier trans- 

form components at the transform frequency v, t\j and tij 
are the positions of the peaks of the two pulse profiles from 
their reference points of Fourier transformation in the time 



domain, and 4>u and <f>2i are the intrinsic phases of the two 
pulse profiles. The cross-spectrum can then be written as, 

CS(u) = F(u)G*(u) = |F(z/)||G(i/)|e _< * csM ; (11) 

where the phase 4>cs (y) is given by 



bcs (v) = 



4>u + 2iri>Atf 



(12) 



with Atf = t2f — tif the fractional time sample bin delay. 

For 4>u{y) = <j>2i(v), i.e. when the pulse profiles at the 
two frequencies have the same shape, the effect of a frac- 
tional bin delay will show up as a linear gradient in the 
phase plot of the CS (see Figure QJ, given by 

A0cs 



At f 



2-kAv 



(13) 



The cross-spectrum can be obtained from the Fourier trans- 
formation of the CCF or from the product of the individual 
Fourier transformations. Of the two we have preferred the 
latter for the AP method, as this helps in the proper prop- 
agation of errors from time domain to frequency domain, as 
explained below. In the SP analysis, however, we have used 
the Fourier transformation of the CCF, with an appropriate 
strategy for computing the errors in the final DM results. 

Let us now look at the estimation of the error in the 
measured delay, which is primarily due to the finite signal to 
noise ratio of the data. For the AP method, the noise in the 
folded profiles, estimated from the off-pulse windows, was 
properly propagated to the CS. For each pulse profile, the 
RMS of phase, c^), and amplitude, o"a(„j, of the Fourier 
transform can be estimated as, 



= &t 



N 



2 (Iml + Rel)' 



(14) 



(15) 



Here N is the number of data points used for fast Fourier 
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transformation (FFT), Im v and Re v are real and imaginary 
parts respectively of the Fourier components at frequency 
bin v, and at is the RMS of the off-pulse noise. The RMS of 
the CS phase, o<$, as M, was computed by adding the noise 
from the two phases in quadrature, 

a lcs(«) - ^iM+%M • ( 16 ) 

In the SP method, the RMS obtained from off-pulse 
windows was properly propagated to estimate the RMS at 
each point of the CCF. The greatest value of this RMS was 
used as a conservative estimate in Equation El to estimate 
the RMS of the CS phase. After this step, the procedure for 
estimating the error in the DM was the same for the AP and 
SP methods. 

The phase gradient, V (4>cs), was computed as the 
slope of the best fitted line, V (bestfit), obtained by the 
least-square method. Thus, 

V (bestfit) x Nfft x T 



At 



360 



(17) 



where Nfft is the number of data points used to compute 
the FFT and T is the time sample. The RMS of At f was 
estimated as 

Cv(bestfit) x Nfft x T 



360 



(18) 



Because the error in At estimation was only due to &At f , 
the error in the final DM value was given by 

0"At, 



At, 



-DM 



(19) 



The above steps were carried out at each epoch to obtain a 
time series of DM values for each pulsar (see Figure [£] for 
example). 



4 RESULTS AND DISCUSSIONS 

The results obtained for the average profile method are 
summarised in Table |5] Here, column 2 gives the catalog: 
DM v alue for each pulsar from lTavlor. Manchester fc Lvnel 
(1993), and the observing frequency bands are given in col- 
umn 3. For each pulsar, we obtained the mean dispersion 
measure over the period of observations, (DM), and the 
quadrature average of C£)M( nojae ) j using: 



(DM) 



2 

&DM, 



2V« 



V^iVep 2 



N, 



where DM, and (7dm, 



(20) 



— . (21) 
are the measured dispersion 



measure and the RMS dispersion at the i th epoch, and N ep is 
the total number of epochs of observations (column 4 of Ta- 
bleEj. The quantity o"DM ( „ oijs) (column 6 of Tabled gives 
the average of the DM error bar estimate from all epochs of 
observations. This quantity gives an estimate of the contri- 
bution to the total RMS fluctuation seen in the time series, 
due to sources of error in the DM estimate. The values for 
c\DM (noisc) for most pulsars are such that the DM estimate 
is accurate to 1 part in 10 4 or better. 

We also estimated the total fluctuation of the DM time 
series, ooM {total) (column 7 of Table[2jl, as 
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Figure 5. Variation of DM with time for pulsars B1642— 03 (up- 
per panel) and B0329+54 (lower panel) observed at frequencies 
610+325 MHz and 610+243 MHz respectively, over the interval 
08 Jan 2001 to 14 May 2002, as a function of day number. The 
continuous line shows the results from average profile analysis, 
and the dotted one from single pulse analysis. The error bars are 
3°'-DM(„ oj5e ) values. 
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Figure 6. DM variation with ScruM , n oite ) error bars for pulsar 
B2217+47 observed at frequencies 610+325 MHz, over the time 
interval 08 Jan 2001 to 14 May 2002 as a function of day number. 
The catalog value of the DM is 43.54 pc/cm 3 . 
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Table 2. DM results from average profile analysis 



Pulsar 


Catalog DM 
(pc/cm 3 ) 


Frequency 
combination (MHz) 




(DM) 
(pc/cm 3 ) 


'BM|„ „ t) 
(pc/cm 3 ) 


a DM (total} 
(pc/cm 3 ) 


ADM/a DM{total) 
(pc/cm 3 ) 


B0329+5 4 


2b. 77b 


243 - 


h 610 


26 


QC 7707(1 

2b. 77870 


0.00003 


n dm no 
0.00103 


+ 2.64 


B0818-13 


40.99 


243 - 


h 325 


32 


40.9222 


0.0013 


0.0043 


- 15.71 


B0823+26 


19.4751 


243 - 


h 325 


29 


19.4545 


0.0004 


0.0016 


- 12.85 




Iz.oo (9 


243 - 


- 325 


29 


iz.8b / 1 


n nfin A 
0.0004 


0.0017 


4- 5.38 


tannin i no 


z.y / Uz 


325 - 


h 610 


•31 


Z.909 ( 


O.OOOo 


O.OOoO 


O 1 
— z. 1 


B1133-t-lb 


4.8471 


325 - 


h 610 


34 


4.8288 


n (Hide 
0.0006 


0.0071 


— 2.57 


B1642-03 


35.665 


325 - 


h 610 


33 


35.75760 


0.00014 


0.00072 


+ 128.20 


B1642-03 


35.665 


243 - 


h 325 


34 


35.72270 


0.00007 


0.00090 


+ 64.00 


B1919+21 


12.4309 


243 - 


h 325 


32 


12.4445 


0.0011 


0.0054 


+ 2.50 


B1929+10 


3.176 


243 - 


h 325 


31 


3.1755 


0.0004 


0.0015 


- 0.31 


B1929+10 


3.176 


325 - 


h 610 


27 


3.1750 


0.0004 


0.0020 


- 0.51 


B2016+28 


14.176 


243 - 


h 320 


29 


14.1611 


0.0007 


0.0025 


- 6.07 


B2016+28 


14.176 


325 - 


h 610 


30 


14.1664 


0.0008 


0.0051 


- 1.90 


B2045-16 


11.51 


243 - 


h 320 


31 


11.5094 


0.0012 


0.0114 


- 0.05 


B2217+47 


43.54 


325 - 


h 610 


31 


43.5196 


0.0007 


0.0061 


- 3.38 



a DM (total) 



(Xgg (M, - (DM)f 



(22) 



In the most general case, this total RMS of the DM fluctu- 
ation is composed of a part due to estimation error on the 
DM (Eaua,tion l21ll and the remaining due to other processes 
likely to play a role in the time variability of DM (a prime 
candidate for which is DM fluctuation due to large scale 
electron density irregularities in the ISM). An estimate of 
the variance due to such processes can be obtained as 



a DM, Isll 



— &DM, 



2 



(23) 



As can be seen in columns 6 and 7 of Table El for almost 
all the pulsars, o-oM {noi3ls) is much smaller than <?r>Mt tota n, 
indicating the presence of substantial DM fluctuations due 
to such sources. We return to this aspect in more detail at 
the end of this section. 



4.1 On the constancy of (DM) estimates 

Keeping in mind the total RMS for each DM estimate 
(o~DM (tatal) ), we can see that the mean DM, (DM), for 
each pulsar is estimated with a fairly good accuracy - 
~ 1 part in 10 3 or better (DM accuracy at each epoch is 
~ 1 part in 10 4 ). It is interesting to compare these mean 
DM values with other estimates in literature. Column 8 
of Table [2] shows the differen ce between our (DM) value 
and the catalog DM value ijTavlor. Manchester fc Lvne I 
Il993h . in units of <JDM, tota n ■ While for most pulsars our 
results agree with the catalog values within 3 (Jnu, tota n , 
there are some pulsars, namely B0818— 13, B0823+26, 
B0834+06, B1642-03 and B2016+28, which show a sig- 
nificant difference. We now discuss these discrepant cases 
in some detail, using for comparison results from ( i) the 
old pulsar catalog of iTavlor Manchester fc Lvne I ll993l l 
(ii) the new pulsar catalog THolTb's^t^rtHl2Tr)4l ~"see also 
www.atnf.csiro.au/research/pulsar/psrcat), and (iii) other 
reports in literature. 

For pulsar B0818-13 we find a (DM) of 40.922 ± 



0.004 pc/cm 3 , which is significantly smaller than the value 
of 40.99±0.03 pc/cm 3 given in the ol d catalog, which comes 
from a very early measurement l Man^h^s^ej^^T^iyjorJ 
Il972h . It is interesting to note that iKuzmin et al. I |1998) 
find an intermediate value of 40.965 pc/cm 3 for this pul- 
sar, from measurements made between 1984 and 1991. Fur- 
thermore, the new pulsar catalog gives a value of 40.938 ± 
0.003 pc/cm 3 iHobhs et al. 12 0041 . w hich is intermediate be- 
tween that of lKuzmin et al. I l|l998h and our result. One in- 
teresting possibility from the above data points is that the 
DM of this pulsar is showing a slow and secular decline with 
time, on time scales of decades. 

For pulsar B0823+26 we again find a mean DM that is 
significantly s maller than the value in t he old catalog (based 



on the work o f lPhillins fc Wolszczan I jl992f)1 . Our result is 
also discrepant from that of Kuzmin et al. I ll998ll . which 



is in good agreement wi th the old catalog v alue. However, 
the new pulsar catalog flHohbs et al. lEooil) cites value of 
19.454 ± 0.004 pc/cm 3 , which is fully consistent with our 
result. 

For pulsar B0834+06 we find a (DM) (12.867 ± 
0.002 pc/cm 3 ) that is somewhat larger than the old cat- 
alog val ue of 12.8579 ± 0.000 2 pc/c m 3 (based on the 
work of IPhillips fc Wolszczan I l|l992h ). For this pulsar, 
IKuzmin et al. I l|l998l ) renort a value of 12.865 pc/cm 3 , 
which agrees quite well with our res ult, whereas the 
new pulsar catalog ijHobbs et al. I 12004 ) cites a value of 
12.889 pc/cm 3 , significantly higher than all the other num- 
bers for this pulsar. 

For B2016+28, our (DM) values (from 2 different pairs 
of frequencies) are consistent with each other, but are signifi- 
cantly smaller than the resu lts cited in the old catalog (based 
on the work of | Craft 11197(1)1 . the new catalo g (based on the 
work of iHobbs et al. I l(2004|l 1. as well as in IKuzmin et al.1 
(1998), all of which are consistent with each other. 



Amongst all our results, the mean value of DM for 
PSR B1642— 03 shows the largest discrepancy with the orig- 
inal catalog v alue of 35.665 ± 0.005 pc/cm 3 (based on very 
early work of iHunt ] l|l97l|l 1. This is true for our DM re- 
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suits from both sets of frequency pairs, though the discrep- 
ancy is more for our results obtained from measurements at 
325 + 610 MHz bands (the difference in DM values from the 
two frequency pairs is discussed separately in the next sub- 
section). We note tha t the DM of 35.73 pc/cm 3 obtained by 
iKuzmin et al. I l)l998l) is equally discrepant from this catalog 
value, and lies in between our two estimates. A similar value 
(35.737 ± 0.003 pc/cm 3 ) is obtained from a multi-frequency 
timing analysis of oyer 30 years of data for this pulsar by 
IShabanova et al~l l)200ll ). The n ew pulsar catalog g ives a 
value of 35.727 ± 0.003 pc/cm 3 ijHobbs et al. 1 120041) . very 
close to the lower of our two results. Clearly, either the orig- 
inal value of the DM reported for this pulsar was erroneously 
estimated, or there has been a significant evolution of the 
DM of this pulsar from the early years of its discovery. 

PSR B1642— 03 is a particularly interesting pulsar, 
in several other respects. There is a significant uncer- 
tainty in the distance estimate to this pulsar. The dis- 
persion measure derived distance is 2.9 kpc (with an un- 
certainty of 50%), whereas the neutral hydrogen measure- 
ments provide a dista nce constraint of 160 pc (lower limit) 
jOraham et al. lll974) . The smaller distance to this pulsar 
is als o supported by a model (by IPrentice and ter Haar I 
1969) that ascribes much of the DM to the presence of the 
Hn region C Oph. alon g the line of sight. Furthermore, 
IShabanova, et al.1 l|200ll) find that this pulsar has a very 
small proper motion and estimate transverse velocities of 
2 and 30 km/s fo r the t wo distance estimates. In addition, 
IShabanova et al.1 feOOlT ) also claim evidence for free preces- 
sion in this pulsar, based on their analysis of the timing 
data. 

Some of the above properties hav e interesting connec - 
tions with the DM results. For example. lOupta et al. I lll994h 
show that the observed scintillation properties of this pul- 
sar are consistent with a line of sight that goes through the 
limb of an Hu region. In such a case, a long term systematic 
variation of the pulsar DM would be expected if there was 
sufficient transverse relative motion between the pulsar and 
the Hu region. However we see no evidence for such a varia- 
tion in our data. On the other and, the observed DM changes 
could be part of a cyclic DM variation on large time scales, 
such as corresponding to the precession period (~ few 1000 
days). For example, results from th e multi-frequency tim- 
ing data of IShabanova et al. I feOOlT ) show timing residuals 
at two different frequencies (~ 0.1 GHz and 0.6—2.3 GHz) 
which have differences between them that vary as a function 
of phase in the precession cycle. The maximum amplitude of 
this difference is ~ 1 ms, implying that the 0.1 GHz pulses 
arrive ~lms later than the pulses at higher frequencies, at 
these phases. One possible explanation of these variations is 
a cyclic change in DM of ~ 2.5 x 10~ 3 pc/cm 3 . This, how- 
ever, is too small compared to the changes and variations 
seen between the different DM values reported above. Thus, 
although there is a lot more information about this pulsar, 
the nature and reason for the observed DM variations does 
not come out clearly. 

From the results in this subsection, it is clear that con- 
stancy of DM estimates (at the level of 1 part in 1000 or 
better) for a pulsar can not be taken for granted. Whether 
these small changes are due to genuine temporal evolution 
of pulsar DMs or due to differences in the estimation tech- 
niques, remains to be established. 



4.2 DM values from different pairs of frequency 
bands 

For two of our pulsars - B1642-03 and B2016+28 - we car- 
ried out the observations at two pairs of frequency bands. 
These data are almost simultaneous in that the observations 
at each epoch were taken within an hour or so of each other, 
and hence can be compared with each other. PSR B1642— 03 
was observed at the frequency pairs of 325 + 610 MHz and 
243 + 325 MHz and we find a significant difference in the 
mean DM values from these two sets of data (Table EJ- 
The value obtained from the higher frequency combination 
(325 + 610 MHz) is higher than that obtained from the lower 
frequency combination (243 + 325 MHz). On the other hand, 
for B2016+28 the DMs obtained from the two frequency 
pairs (325 + 610 MHz and 243 + 320 MHz) are the same 
within errors, as determined by &DM, total y 

Though it is generally thought that the DM value for a 
pulsar is independent of the frequency of measurement, there 
have been reports in literature about differences in pulsar 
DMs that have bee n estimated from di ff erent parts of t he 
radio spectrum fe.g. lShitov et al. I l|l98sf ): lHankins1 Jl99 Jl ). 
In most of these results, the evidence is for an excess delay in 
the arrival of the pulses at low frequencies, when attempting 
to align them with a DM value computed from the higher 
frequencies. However, our results for PSR B1642— 03, albeit 
for a relatively narrow range of radio frequencies, show an 
opposite trend in that the DM value is larger for the higher 
frequency pair (325 + 610 MHz). 

There are different possible explanations for frequency 
dependent DM variations. For example, an evolution in the 
shape of the profile with frequency can play a role in chang- 
ing the inferred alignment between the profiles at two differ- 
ent frequencies. This should play a more significant role for 
pulsars with complex, multi-component profiles, but should 
be relatively insignificant for pulsars with simple profiles 
(such as pulsars B1642— 03 and B2016+28 in our sample). 
Another interesting possibility is an extra time delay be- 
tween emission received at two frequencies due to differ- 
ent heights of emissi on of these frequencies in the pulsar 
magneto-sphere fe.g. iKardashev et al. lll982T ). an idea that 
has not received much attention in the past. These aspects 
will be examined in greater detail in a separate, forth-coming 
paper. 



4.3 DM differences from average profiles and 
single pulses 

As described in section the DM estimates were obtained 
from two independent methods: measurement of delays be- 
tween the average profiles (the AP method; results reported 
in Table [21 and measurement of delays between single pulse 
trains (the SP method). We find, in general, that the DM re- 
sults for a pulsar depend on the method of analysis. For some 
pulsars, this difference is negligible, e.g. PSR B1642— 03 
(see Figure 0. For others it is significant: PSR B0329+54 
is one such example in our study (see lower panel of Fig- 
ure - the (DM) value obtained from the SP analysis is 
26.7751 ± 0.0007 pc/cm 3 , which is significantly lower than 
the catalog value, which in turn is lower than the (DM) 
value from the AP analysis. 

It is worth noting that dispersion measure values es- 
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timated from alignment of average profiles and those from 
cross-correlation of single pulse emission fe atures have been 
reported to be different in the past also. lHankinsI ilflfllT l 
found that DM values from average profile measurements 
are significantly larger than those obtained from cross- 
correlation of pulsar micro-structure, for P SR B0950+08 
and PSR B1133+16. IStinebring et al. I jl992h have also in- 
vestigated results for PSR B1133+16 over a ten year period, 
obtained using different techniques, and found significant 
variations in the DM values. 

Further, as described in section El these two methods 
actually measure slightly different quantities. Thus, the dif- 
ference between average profile and single pulse analysis re- 
sults that we find is not so surprising. A detailed description 
of these results and an investigation into the possible causes 
and implications of the same will be taken up in another 
forthcoming paper. 

4.4 Slow fluctuations of pulsar DMs 

As mentioned at the beginning of this section, there is evi- 
dence for substantial temporal fluctuations in DM values for 
most of the pulsars. A large part of this is likely to be due 
to the ISM. A detailed study of this aspect will be taken up 
separately in another forthcoming paper. Here, we briefly 
comment on the variations observed, comparing them with 
earlier published results. 

Variations in pulsar DM, by definition, can arise due 
to either spatial and temporal changes in the electron den- 
sity along the line of sight, or change in the distance to 
the pulsar, or both. Electron density changes along the line 
of sight to the pulsar can be in the form of fluctuations 
resulting in DM fluctuations; alternatively, there can be a 
monotonic increase (or decrease) in DM due to the pulsar 
sampling a gradient of the electron density. Most of our ob- 
served DM fluctuations (except PSR B2217+47) show fluc- 
tuations over a constant mean DM, indicating that the ob- 
served changes are due to electron density fluctuations in 
the ISM. In the case of temporal change of distance to the 
pulsar with respect to the observer, the effect would man- 
ifest only as monotonic increase or decrease in pulsar DM. 
In our sample, PSR B2217+47 shows a monotonic increase 
of its DM (FigureEJ- However, the amplitude of this change 
(~ 0.02 pc/cm 3 /year) is such that it would require a very 
large radial velocity (~ 10 6 km/s) through a normal density 
ISM (~ 0.02 /cm 3 ), or a very high density ISM (~ 200 /cm 3 ) 
for normal pulsar velocities (~ 100 km/s). It is likely that 
the cause for this change is due to the pulsar sampling an 
electron density gradient in the ISM, rather than due to 
radial motion of the pulsar. 

Such an electron density gradient can be produced by 
the line of sight to the pulsar crossing through a blob of 
enhanced plasma density. Taking the electron density en- 
hancement of An e pc/cm 3 in a wedge of thickness L pc, 
the observed change in DM is ADM = An e .L . The pul- 
sar's transverse displacement X (of ~ 3 x 10~ 4 pc), sam- 
ples this electron density gradient in one year. Assuming 
the wedge to be part of a spherical blob of radius L pc 
(and X ~ L), we can estimate the electron density gradi- 
ent to be ~ 2 x 10 5 /cm 3 /pc or 1 /cm 3 /AU. This value is 
a lower limit - if the cloud is closer to the observer, the 
electron density gradient could be even higher. Evidence 



for such AU-size clouds of enhanced electron density in the 
ISM also co mes from scint i llatio n observations of pulsars. 
For example, iBhat et. al. I lll999f) find evidence for clouds 
with length scales of ~ 10 AU and electron density contrast 
~ a few electrons/cm 3 . Our results are similar, though a bit 
on the higher side. 

Long term, slow DM variations, on time scales of weeks 
to mo nths, have been st udied in the past bv [ B acker e t al. I 
<199fll (3 pulsars) and IPhillins fc Wolszczan I il99lh (6 
pulsars). Whereas iBacker et al. I il993h report total DM 
fluctuations ~ 0.02 p c/cm 3 over 1 — 2 year periods, 
|p"h illins fc Wolszczanl l|l99lfl report typical variations ~ 
0.002 pc/cm 3 (and smaller) over similar time intervals. 
Our results show aDM {total) ~ 0.001 to 0.007 pc/cm 3 
for most cases, implying total fluctuations ~ 0.005 to 
0. 03 pc/cm 3 . These are typicall y larger than those reported 
bv lPhill ips Wolszczan | | l991fl . but comparable to the re- 
sults of Backer et al. I jl993h . 

4.5 Summary 

We have presented a new experiment for accurate measure- 
ment of pulsar DMs using the GMRT in a simultaneous, 
multi-frequency sub-array mode. We have shown that single 
epoch DM estimates using this technique can achieve an ac- 
curacy of 1 part in 10 4 or better. With improved sensitivity 
performance of the GMRT and faster sampling that is now 
available, this accuracy can be improved in future experi- 
ments and the technique can be extended to a larger set of 
pulsars. From the large number of epochs of DM measure- 
ments for each of the 12 pulsars in our sample, we are able to 
obtain fairly accurate estimates for the mean DM for most of 
them. A detailed comparison of DM values in the literature 
with our mean DM values highlights the lack of consistency 
(at the level of ~ 1 part in 1000) in the different DM esti- 
mates, the reason for which remains to be understood. We 
have also briefly highlighted some of the other results from 
our data - such as DM estimates from different frequency 
combinations, differences in average profile and single pulse 
DM values, and slow fluctuations of pulsar DMs with time 
(which are most likely to be due to ISM effects) - these will 
be the subject of follow-up papers. 
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